American Journal of Business Education - December 2009 


Volume 2 , Number 9 


A MATLAB-Aided Method For Teaching 
Calculus-Based Business Mathematics 

Jiajuan Liang, University of New Haven, USA 
William S. Y. Pan, University of New Haven, USA 


ABSTRACT 

MATLAB is a powerful package for numerical computation. MATLAB contains a rich pool of 
mathematical functions and provides flexible plotting functions for illustrating mathematical 
solutions. The course of calculus-based business mathematics consists of two major topics: 1) 
derivative and its applications in business; and 2) integration and its applications in business. 
Both topics involve various mathematical functions that may scare many business students. 
Fortunately, MATLAB provides an easy solution to handling the derivative and integration of 
various mathematical functions. With the help of MATLAB, all elementary formulas for derivative 
and integration in the textbooks of calculus-based business mathematics become less important to 
those students who have difficulty in learning and using mathematical formulas to find solutions to 
quantitative business problems. With the wide availability of various business data bases 
including the MATLAB package in the Samuel S. Bergami, Jr. Learning Center for Finance and 
Technology at the University of New Haven, in this paper, we illustrate how to use MATLAB to 
help students’ learning of calculus-based business mathematics and apply the calculus principles 
to find solutions to real quantitative business problems. Our experience shows that MATLAB 
greatly saves students’ efforts in learning difficult and pure calculus formulas and it is easy for 
business students to acquire the basic internal functions in using MATLAB. 


1. INTRODUCTION 



alculus-based business mathematics is a difficult course to many business students due to their poor 
algebraic skills. It is also challenging for instructors to teach this course within a limited time schedule, 
usually one-semester for covering derivative and integration topics, and a small part of multivariate 
calculus and applications in business. Without resorting to modern technology, traditional teaching of this course 
may be just presenting the formulas to students without proof, illustrating by examples and then having students do 
simple exercises by hand. This old-style teaching method may receive desirable effects if all students have the 
necessary preliminary algebra knowledge and algebra-operation skills. Many business students, however, show 
strong dislike in mathematical formulas. Some are even math-phobia when doing algebraic operations with 
relatively complicated mathematical functions by hand. It is a common situation that students in the same class have 
various levels of algebraic knowledge. When having students solve derivative or integration problems, after the 
instructor’s illustrations of using the formulas, some students may be able to follow the same idea in the illustrated 
examples, but others have no idea in how to start with —they simply couldn’t get the idea from the instructor’s 
illustration. When facing with this situation, the instructor may either repeat the illustrations a few more times or 
make up the missing algebraic skills for those students in order to move on to the next topics. This may not be a 
problem if the teaching time is not limited. However, business mathematics is usually a one-semester course to 
business students and regular topics in derivative and integration and their applications in business are supposed to 
be covered. This calls for the resort of modern technology to help students’ learning. Widely available and 
affordable computer programs make it possible for business students to learn calculus-based business mathematics 
with substantially reduced efforts, especially for those students with poor or weak math-operation skills. 


Among the widely available symbolic computer programs such as MATHEMATICA, MATLAB and 
others, MATLAB may be one of the easiest one for business students to learn in short time. The Symbolic Math 
Toolbox in MATLAB provides almost all necessary math-operation techniques for students to the course of business 
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mathematics. It is a quite direct approach to learning MATLAB commands for math operations in calculus, 
optimization and integration. In this paper, we will summarize our experience in using MATLAB to teach business 
mathematics for undergraduate students in the College of Business at the University of New Haven. Section 2 
presents the step-by-step MATLAB-aided teaching topics in our one-semester course of business mathematics. 
Some concluding remarks are given in the last section. 

2. MATLAB-AIDED TEACHING TOPICS 


2.1. Introduction to limit 


The major topics in calculus are based on the concept of limit. Most students learned this concept in high 
school but some of them forgot most of the mathematical operation in obtaining the limit. To proceed to the concept 
of instantaneous change rate or derivative of a function, it is necessary to review the concept of limit and show 
students how to obtain the limit of a given function. Because the topic of limit is transitory in teaching calculus, we 
just review the concept and show the students how to use MATLAB to obtain the limit of a given function. The 
following example covers the limits of frequently used functions such as polynomial functions, rational functions, 
root functions, composite power functions, exponential and logarithmic functions. 


Example 2.1. Computing the limits. 

(1) limh 3 — 2x — \) 

Jt->2 V ’ 


(3) lim 


*->i 


A 2 -i 

x-\ 


+ V x 2 + 3 


(5) lim4o(l-e“° 12f ) 

r ->10 v ' 


(2) lim 

x —>3 


x + 3 


- + 


x-3 
x 2 -9 


\ 

) 


(4) lim(l + 2x) llx 

x->0 

(6) lim[ln(2x - 2) - ln(x -1)] 


By using the “syms” and “limit” commands in MATLAB, the output can be obtained: 


( 1 ) 



This is a Classroom License for instructional use only. 
Research and commercial use is prohibited. 

To get started, select MATLAB Help or Demos from the Help menu. 

» syms x 

» limit('x A 3-2*x-1',x,2) 


( 2 ) 


MATLAB 7.3.0 (R2006b) 


File Edit Debug Desktop Window Help 

□ Q 1 | f> current Directory: [ C:\Program Files\MATLAB\R2006b\work 

Shortcuts 0 How to Add 0 What's New 


This is a Classroom License for instructional use only. 
Research and commercial use is prohibited. 

To get started, select MATLAB Help or Demos from the Help menu. 

» syms x 

» limit( , x/(x+3)+(x-3)/(x A 2-9)',x, 3) 


ans = 


2/3 


3 
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(3) 



This is a Classroom License for instructional use only. 
Research and commercial use is prohibited. 

T o get started, select MAT LAB Help or Demos from the Help menu. 

» syms x 

» limit('(x A 2- 1 )/(x-1 )+sqrt(x A 2+3) , ,x, 1) 
ans = 

4 

(5) 



This is a Classroom License for instructional use only. 
Research and commercial use is prohibited. 

To get started, select MATLAB Help or Demos from the Help menu. 

» syms t 

» limit('40*(1-exp(-0.12*t))',t,1 0) 
ans = 

27.952231523511916134200895716671 


(4) 



This is a Classroom License for instructional use only. 
Research and commercial use is prohibited. 

To get started, select MATLAB Help or Demos from the Help menu. 

» syms x 

» limited +2*x) A (1/x)',x,0) 

ans = 
exp (2) 


» exp (2) 

ans = 

7.3891 


( 6 ) 



This is a Classroom License for instructional use only. 
Research and commercial use is prohibited. 

To get started, select MATLAB Help or Demos from the Help menu. 

» syms x 

» limiteiog(2*x-2)-log(x-1)',x,2) 

ans = 

l°g(2) 

»log(2) 
ans = 


0.6931 


2.2. Graphing and its application 

Graphing is an intuitive and the easiest way for students to understand the numerical relationship between 
two variables or to compare the increasing and decreasing rates of two or more functions simultaneously. The 
instantaneous change rate of derivative of a function can be easily interpreted by graphing. 

Example 2.2. Using the graphing technique in identifying the optimal business strategy. 

(1) A company manufactures and sells v transistors per week. The weekly cost and price-demand equation are 
given by 
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Cost function: C(x) = 90,000 + 30x, 0 < X < 9,000, 

x 

And the price function: p = 300 — — . Identify the approximate optimal demand x to maximize the profit. 
The profit function is obtained: 

P(x ) = Revenue - Cost = 
xp - C(x) = x 


f 

300- 


v 


30 


(90,000 + 30x) =-+ 27Ox - 90,000, 0 < x < 9,000. 

30 


By using the “ezplot” command in MATLAB, the output is obtained. It can be easily identified that the approximate 
optimal demand x to maximize the profit is close to 4,000 transistors. 



File Edit Debug DeskJtop Window Help 


□ & 

I A % (Hi 

* d □ 

f 

Current Directory: 

C:\Program Files\MATLAB\R2006b\work v [27] [t] 


Shortcuts [7] How to Add 0 What's New 

This is a Classroom License for ir 
Research and commercial use is 

To get started, select MATLAB Help or 

» ezplot('-x A 2/30+270*x-9000 l , [0,9000]) 

» 


Figure 1 


File Edit View Insert Tools DeskJtop Window Help ^ 

□ » aa 7 I 7T© "W □’b * 


-x 2 /30+270 x-9000 



Figure 1. Plot of the profit function versus the demand level 


(2) A mail-order company specializing in computer equipment estimated (based on past data) the price-demand 
equation 

p = exp (5.4474 - 0.0008x), 200 < x < 900 

Where x is the weekly demand for x modems with unit price $100. Identify the approximate optimal demand x to 
maximize the profit. 

The profit function is obtained: 

P(x) = Revenue - Cost = xexp(5.4474 - 0.0008x) - lOOx, 200 < x < 900. 

By using the “ezplot” command in MATLAB, the output is obtained. It can be easily identified that the approximate 
optimal demand x to maximize the profit is close to 450 modems. 
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This is a Classroom License for instructional u 
Research and commercial use is prohibited. 

To get started, select MATLAB Help or Demos from 

» ezplot('x*exp(5. 4474-0.0008*x)-1 00*x\[200,900]) 

» 


Figure 1 




File Edit View Insert Tools Desktop Window Help 


□ & y 


^c?® « ns 


x10 4 x exp(5.4474-0.0008 x)-100 x 



Figure 2. Plot of the profit function versus the demand level 

(3) A department store estimated the (based on past data) the price-demand equation for selling x cream puffs 
as 

p = 6.2242 - 0.5238 ln(x), 5,000 < jc < 10,000. 

Assume the unit cost for making a cream puff is $1. Identify the approximate optimal demand x to maximize the 
profit. 

The profit function is obtained: 

P(x) = Revenue - Cost = x(6.2242 — 0.5238ln(x)) — x, 5,000 < x < 10,000. 

By using the “ezplot” command in MATLAB, the output is obtained. It can be easily identified that the approximate 
optimal demand x to maximize the profit is close to 8,000 cream puffs. 



This is a Classroom License for instructional use only. 
Research and commercial use is prohibited. 


To get started, select MATLAB Help or Demos frorr 


» ezplot('x*(6.2242-0.5238*log(x))-x',[5000,1 0000]) 

» 


File Edit View Insert Tools DeskJtop Window Help 


□ y 


^ ^ n ® 


x (B.2242-0.5238 log(x))-x 



3800 

5000 6000 7000 8000 9000 10000 


Figure 3. Plot of the profit function versus the demand level 
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Example 2.3. Graphing two or more functions in a single window. Use the graphing technique to identify the 
demand interval in which profit can be made. When the revenue is greater than the cost, profit can be made. 
Therefore, it is necessary to plot the revenue function and cost function in a single window for comparison between 
the two functions. 


(1) In (1) of Example 2.2, the two functions are: 


Revenue function : R(x) = x 


f 


300-- 


30 


Cost function : C(x ) = 90,000 + 30v, 0 < v < 9,000. 


The “plot” command in MATLAB is used to plot two or more functions in a single window. This command requires 
generating the plotting points first before plotting. The MATLAB steps are given as in the following output. 

Step 1. Generate the plotting values for v. Under the MATLAB command sign, type “x=[0:10:9000];”, where “0” is 
the starting point of the plot and “9000” the ending point, “10” is the step length which is chosen by convenience; 

Step 2. Generate the plotting values for C(x) by the given cost function by typing “C=90000+30*x;” under the 
MATLAB command line, here the semi-colon “;” is typed to avoid any output from the command; 

Step 3. Generate the plotting values for R(x) by the given revenue function by typing “R=x.*(300-x/30);” under 
the MATLAB command line, where the “dot product” “.*” is used for product of vectors because x=[0:10:9000] is a 
vector of 1x900; 


Step 4. After generating the plotting values, it is ready to plot by using the “plot” command. In the “plot” command, 
“x,R,’-‘” means that the plot for R(x) versus v is a “real line”, “x,C,’—‘” means that the plot for C(x) versus v is 
a “dashed line”. MATLAB provides several options for the types of line; 

Step 5. Use the “gtext” command to add a label to the curves for R(x) and C(x) , respectively. This step is just 
optional. Type “gtext(‘... ’)” under the MATLAB command line for whatever text “...” you want to label. 

Step 6. Find the break-even points (zero profit): the revenue is equal to the cost so that zero profit is obtained by 
using the “solve” command in MATLAB to solve the equation “ R(x) = C(x) ”, which gives two solutions 

X, =4050-150 x609 1/2 «348 and x 2 = 4050 +150 x 609 1/2 * 7752. 


Conclusion: when the demand level is 348 < X < 7752, it makes positive profit. Otherwise, it results in negative 
profit. The break-even demand levels are x=348 and x=7752 , where zero profit is obtained. 
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This is a Classroom License for instructional use only. 
Research and commercial use is prohibited. 

To get started, select MATLAB Help or Demos from the Help menu. 


»x=[0:1 0:9000]; 

» C=90000+30*x; 

» R=x.*(300-x/30); 

» plot(x, R j'-'jXjC,'— ') 

» gtext('R(x)') 

» gtextCC(x)') 

» solve('x*(300-x/30)=90000+30*x') 
ans = 

4050-1 50*609 A (1/2) 

4050+1 50*609 A (1 /2) 



Figure 4. Plot of the revenue function and cost function versus the demand level 


(2) In (2) of Example 2.2, the two functions are: 

Revenue function : R{x) = xp = x exp (5.4474 - 0.0008x), 

Cost function : C(x) = lOOx, 200 < x < 900. 

By following the same steps as above, the MATLAB output is given as follows. It is clear tat the break-even 
solutions are Xj = 0 and x 2 « 1053. 

Conclusion: when the demand level is 0 < X < 1053, it makes positive profit. Otherwise, it results in negative 
profit. The break-even demand levels are x=0 and x=1053 , where zero profit is obtained. 
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This is a Classroom License for instructional use only. 
Research and commercial use is prohibited. 


To get started, select MAT LAB Help or Demos from the Help menu. 


» x=[200:10:1600]; 

» C=100*x; 

» R=x.*exp(5.4474-0.0008*x); 

» plottx.R.'-'.x.C,'—') 

» gtext('R(x)') 

» gtext('C(x)') 

» solve('x* exp(5.4474-0.0008*x)=100*x') 


ans = 


0 . 

1052.7872675148857899550213632891 



Figure 5. Plot of the revenue function and cost function versus the demand level 


(3) In (3) of Example 2.2, the two functions are: 

Revenue function : R(x) = xp = v(6.2242 - 0.5238ln(v)), 

Cost function : C(x) = x , 5,000 < v < 10,000. 

By following the same steps as above, the MATLAB output is given as follows. It is clear tat the break-even 
solution is X « 21,454 . 

Conclusion: when the demand level is X < 21,454 , it makes positive profit. Otherwise, it results in negative profit. 
The break-even demand level is x=21,454 , where zero profit is obtained. 
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This is a Classroom License for instructional use only. 
Research and commercial use is prohibited. 

To get started, select MAT LAB Help or Demos from the Help menu. 


» x=[5000:10:30000]; 

» C=x; 

» R=x.*(6.2242-0.5238*log(x)); 

» plot(x,R,'-\x,C,'~') 

» gtext('R(x)') 

» gtext('C(x)') 

» solve('x* (6.2242-0.5238*log(x))=x') 
ans = 

21453.735672905700785456602748 



Figure 6. Plot of the revenue function and cost function versus the demand level 


2.3. Derivative and optimization 

2.3.1. The first derivative and its application in business 

Derivative is the major topic in teaching calculus-based business mathematics. The major application of 
calculus in business is to find the optimal solution to a quantitative business problem. This is solved by optimization 
process based on derivative. Therefore, derivative is the key to solve many quantitative business problems related to 
optimization. In widely available textbooks on business mathematics, a common way to handle the derivative topic 
is to illustrate by various examples how to use the formulas to obtain the derivative of a given function. Without 
resorting to modern computer programs, this method may work for students with sufficient algebra-operation skills. 
Business students are usually more interested in finding solutions to business problems without preparing 
themselves in basic mathematical theory. Resorting to modern computer programs such as MATLAB with less 
emphasis on manual operation makes it possible for students to realize this purpose. By using MATLAB, most 
students realized that it is a lot of fun instead of burden to study business mathematics with the help of MATLAB. 
We illustrate the use of MATLAB in teaching derivative topics through the following example. 

Example 2.4. Price-supply and price-demand problems. 

(1) The number of v stereo speakers a retail chain is willing to sell per week at a price of $p is given by 


x = 80 ^/p + 25 - 40, 20 < p < 100. 
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Find the derivative x'(p) = - and x' (75) , and interpret the result. 

dp 

(2) The number of y stereo speakers people are willing to buy per week from a retail chain at a price of $p is 
given by 

y = 1000- 60/p + 25, 20 < p < 100. 


Find the derivative y' (p) = — ^ and (75) , and interpret the result. 

dp 


(3) What is the common number of stereo speakers that retail chain is willing to sell and people are willing to 
buy per week? 

The solutions are given by MATLAB output in Figure 7 as follows. 


(1) Derivative computation 



This is a Classroom License for instructional use only. 
Research and commercial use is prohibited. 

To get started, select MATLAB Help or Demos from the Help menu. 

» diff( , 80*sqrt(p+25)-40') 


40/(p+25) A (1/2) 


(2) Derivative computation 


MATLAB 7.3.0 (R2006b) 


File Edit Debug Desktop Window Help 


Da; HTh m ^ » d o 

^ Current Directory: C:\Program Files\MATLAB\R2006b\work s 

Shortcuts 0 How to Add 0 What's New 



This is a Classroom License for instructional use only. 
Research and commercial use is prohibited. 

To get started, select MATLAB Help or Demos from the Help menu. 

» diff('1 000-60*sqrt(p+25)') 


-30/(p+25) A (1/2) 


» p=75; 

» 40/(p+25) A (1/2) 
ans = 

4 


x\p) = 40/(/? + 25) 1/2 , x'(75) = 4. 


Figure 7 a 


» p=75; 

»-30/(p+25) A (1/2) 
ans = 

-3 

y’(p) = -30/(p + 25) 1/2 , y’C 75) =-3. 

Figure 7b 
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This is a Classroom License for instructional use only. 
Research and commercial use is prohibited. 

To get started, select MATLAB Help or Demos from the Help menu. 


» solve( , 4*75+b=80*sqrt(75+25)-40') 


ans = 


460 


» p=[20:1:100]; 

» x=80*sqrt(p+25)-40; 

» y1=4*p4460; 

» plot(p,x, , - , ,p,y1, , -- , ) 

» gtext('x= 8 0 * s q rt(p+25)-40 ') 

» gtext('tangent line yl =4p4460‘) 

» gtext('at p=75, slope=dx/dp(75)=4') 

Figure 7c Interpretation of x (75) = 4 


fc-Jigure L 00® 



MATLAB 7.3.0 (R2006b) 


File Edit Debug Desktop Window Help 


D \2 f j n* j ^ ^ ^ Current Directory: C:\ProgramFiles\MATLAB\R2O06b\work v [7 

3 * 

Shortcuts 0 How to Add [7] What's New 



This is a Classroom License for instructional use only. 
Research and commercial use is prohibited. 


To get started, select MATLAB Help or Demos from the Help menu. 


» solve('-3*75+b=1 000-60*sqrt(75+25)') 


ans = 


625 


» p=[20:1:100]; 

» y=1000-60*sqrt(p+25); 

» y2=-3*p+625; 

» plot(p,y, , - , ,p,y2, , - , ) 

» gtext('y=1 00Q-60*sqrt(p+25)') 

» gtext('tangent line y2=-3p+625‘) 

» gtext('at p=75, slope=dy/dp(75)=-3‘) 



Figure 7d Interpretation of y'(75) = — 3 
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(3) the common number of stereo speakers that retail chain is willing to sell and people are willing to buy per week 
is the intersection point of the curve for A = 80 ^/p + 25 — 40 and the curve for y = 1000 - 60-7 P + 25 • That 


is, 

80-7 /? + 25 - 40 = 1000 - 607 p + 25 


Which gives the solution x = y « 554 stereo speakers at the common price p — 1479/49 « $30.18 per unit. 



This is a. Classroom License for instructional use only. 
Research and commercial use is prohibited. 

To get started, select MATLAB Help or Demos from the Help menu. 

» solve('80*sqrt(p+25)-40=1 000-60*sqrt(p+25) , > 

ans = 

1479/49 


» p= 1 479/49; 

» 80*sqrt(p+25)-40 

ans = 

554.2857 

Figure 7e MATLAB Computation of the common price and the common number of stereo speakers that retail chain is 

willing to sell and people are willing to buy per week 


Example 2.5. (Example 2.2 continues) Maximizing profit, marginal analysis in business and economics. Find the 
marginal profit and the best demand to maximize the profit for the profit functions in Example 2.2. 


( 1 ) 


The profit function is: 


P(x) = x 


f 

300 - — 

V 30y 


-(90,000 + 30x), 0<x< 9,000. 


Marginal profit is the derivative P } (a) , and the best demand satisfies P ? (a) = 0. The results are given in the 


MATLAB output in Figure 8, which gives P'(a) = 270 — —, and the best demand x=4050 that satisfies 
P' (a) = 0. Figure 1 in Example 2.2 shows that x=4050 really leads to the maximum profit. 
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This is a Classroom License for instructional use only. 
Research and commercial use is prohibited. 

To get started, select MATLAB Help or Demos from the Help menu. 

» diff( , x*(300-x/30)-(90000+30*x)') 

ans = 

270-1/1 5*x 


» solve('270-1 /I 5*x=0‘) 
ans = 

4050 

Figure 8. MATLAB computation of marginal profit and the best demand to maximize the profit 

(2) The profit function is: 

P(x ) = x exp (5.4474 - 0.0008x) — lOOx, 200 < x < 900. 

The results are given in the MATLAB output in Figure 9, which gives the marginal profit 
P'(x) = exp(5.4474 — 0.0008x) -0.0008xexp(5.4474 - 0.0008x) -100, and the best demand x-467 
that satisfies P' (x) = 0 . Figure 2 in Example 2.2 shows that x=467 really lead to the maximum profit. 



This is a Classroom License for instructional use only. 
Research and commercial use is prohibited. 

To get started, select MATLAB Help or Demos from the Help menu. 

» diff('x*exp(5. 4474-0.0008*x)-1 00*x') 

ans = 

exp (5.4474-.8e-3*x)-.8e-3*x*exp(5.4474- .8e-3*x)-1 OO 


» solve('exp(5 .4474- 8e-3*x)-.8e-3*x*exp(5.4474- ,8e-3*x)-100=0') 
ans = 

467.41573463086774660112793277537 
Figure 9. MATLAB computation of marginal profit and the best demand to maximize the profit 
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P(x) = x(6.2242 - 0.5238 ln(x)) - x, 5,000 < x < 10,000. 

The results are given in the MATLAB output in Figure 9, which gives the marginal profit 
P\x) = 4.7004 —0.5238 ln(x), and the best demand x=7892 that satisfies P\x ) = 0. Figure 3 in Example 
2.2 shows that x=7892 really leads to the maximum profit. 



This is a. Classroom License for instructional use only. 
Research and commercial use is prohibited. 

To get started, select MATLAB Help or Demos from the Help menu. 

diffOc* (6.2242-0.5238*log(x))-x') 

ans = 

4.7004- .5238*log(x) 


solve ('4 .7004-.5238*log(x) = 0') 
ans = 

7892.3882903883863013266956575602 


Figure 10. MATLAB computation of marginal profit and the best demand to maximize the profit 
2.3.2. The second derivative and its application 

The major application of the second derivative in business is to confirm whether the optimal solution to the 
problem by the first derivative method leads to maximum or minimum. This requires checking the sign of the 
second derivative at the optimal point(s). A negative value of the second derivative at an optimal point confirms that 
this optimal point is a local maxima; a positive value of the second derivative at an optimal point confirms that this 
optimal point is a local minima. Usually speaking, a local optimal solution is also the absolute optimal solution to a 
business problem. 

Example 2.6. Example 2.5 continues. 

(1) In Example 2.5, the optimal solution to the profit function is x=4050. The computation of the second 
derivative of the profit function at this point is given in Figure 11, which is a negative constant. This 
confirms that x=4050 is a local maxima. This implies that the optimal demand level is x=4050. 
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This is a Classroom License for instructional use only. 
Research and commercial use is prohibited. 

To get started, select MATLAB Help or Demos from the Help menu. 

» diff0**(300-x/30)-(90000+30*x) , ,2) 

ans = 

-1/15 


Figure 11. Computation of the second derivative of the profit function 

(2) In Example 2.5, the optimal solution to the profit function is x-467. The computation of the second 
derivative of the profit function at this point is given in Figure 12. 

P"( X) = -0.0016 exp (5.4474 - 0.0008x) + 0.64 x 10^ x exp (5.4474 - 0.0008x) 

And P M (467) = —0.2079 < 0. This confirms that x=467 is a local maxima. This implies that the optimal 
demand level is x-467. 



This is a Classroom License for instructional use only. 
Research and commercial use is prohibited. 

To get started, select MATLAB Help or Demos from the Help menu. 

» diff('x*exp(5.4474-0.0008*x)-1 OCTx',2) 

ans = 

-. 1 6e-2*exp(5.4474- .8e-3*x)+.64e-6*x*exp(5.4474- .8e-3*x) 


» x=467; 

>> -.1 6e-2*exp(5.4474-.8e-3*x)+.64e-6*x*exp(5.4474- .8e-3*x) 
ans = 

-0.2079 


Figure 12. Computation of the second derivative of the profit function 


(3) 


In Example 2.5, the optimal solution to the profit function is x=7892. The computation of the second 
derivative of the profit function at this point is given in Figure 12. P"(x) = — 0.5238 /x and 
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P"(7 892) = — 6.6371 xlO -5 < 0. This confirms that x-7892 is a local maxima. This implies that the 
optimal demand level is x=7892. 



This is a Classroom License for instructional use only. 
Research and commercial use is prohibited. 

To get started, select MATLAB Help or Demos from the Help menu. 

» diff( , x*(6.2242-0.5238*log(x))-x , ,2) 

ans = 

-5238/x 


x=7892; 

-.5238/x 

ans = 

-6.6371 e-005 

Figure 13. Computation of the second derivative of the profit function 


2.4. Integration and its applications in business 

Integration may be the most challenging topic in teaching calculus-based business mathematics for both 
instructors and students. Traditional teaching of this topic is usually by presenting the formulas and illustrating the 
application of the formulas through various examples. This again requires students’ even stronger algebra-operation 
skills. It is just these skills that make it the most difficult topic for students to learn integration. Integration consists 
of indefinite integrals and definite integrals. While the idea for both indefinite and definite integrals is usually easy 
for business students to understand, manual computation of the indefinite and definite integrals is completely 
another story for many students. Based on our experience, MATLAB makes it possible for business students to learn 
the integration topic and apply it to business problems without much difficulty. 

Example 2.7. Revenue analysis. The weekly marginal revenue from the sales of v pairs of tennis shoes is given by 

/?'(x) = 40-0.02*+ 222, /?(0)=0, 

x + 1 


Where R(x ) is the revenue in dollars. Find the revenue function and the revenue from the sale of 1,000 pairs of 
tennis shoes. 

The MATLAB solution is given in Figure 14, where the first command gives the indefinite integral 


R(x) = J( 40 -0.02* 


+ 


200 
x + 1 


\ 

dx = 40x - O.Obc 2 + 2001n(x +1) + C 


J 
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(C does not appear in the MATLAB output). The second command is for identifying the constant C from the initial 
condition /?(0) = 0 , which gives C=0. So the final solution is 


R(x) = J[40 - 0.02x + 


200 N 
X + 1 


dx = 40x - O.Olx 2 + 2001n(x + 1) with /?(0) = 0. 


It is easy to obtain 

R(l 000) = 40(1000) -0.01(1000) 2 + 200ln(l000 + 1) = $31,382. 



This is a Classroom License for instructional use only. 

Research and commercial use is prohibited. 

To get started, select MATLAB Help or Demos from the Help menu. 

» int( l 40-0.02*x+200/(x+1)') 

ans = 

40.*x-.1 OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOe-1 *x A 2+200.*log(x+1 .) 

» solveO40.*0-.1 Oe-1 *0 A 2+200.*log(0+1 .)+C=0') 
ans = 

0 


Figure 14. Computation of indefinite integral with a given initial condition 

Example 2.8. Marketing analysis. An automobile company is ready to introduce a new line of cars with a national 
sales campaign. After test marketing the line in a carefully selected city, the marketing research department 
estimates that sales (in millions of dollars) will increase at the rate of 

S'(t) = 10-lOe” 01 ', 0 < ? < 24, 

t months after the national campaign has started. 

(1) What will be the total sales, S(t) , t months after the beginning of the national campaign if assuming zero sales 
at the beginning of the campaign? 

(2) What are the estimated total sales for the first 12 months of the campaign? 

(3) When will the estimated total sales reach $100 million? 
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The same steps as in Example 2.7 are followed and the MATLAB output for (1) is given in Figure 15a. It shows that 


Sit) = 10* + 100exp(-0.1/) -100. 



This is a Classroom License for instructional use only. 
Research and commercial use is prohibited. 

To get started, select MATLAB Help or Demos from the Help menu. 

int(‘ 1 0-1 0*exp(-0.1 *t)') 

ans = 

1 Q.*t+1 0O.*exp(-.1 OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOT) 


solve('1 0 *0 + 1 00.*exp(-.1 *0)+0 = 0 , > 
ans = 


-1 OO. 


Figure 15a Computation of the indefinite integral with a given initial condition 

Results for (2) and (3) are given in Figure 15b, where for (3), it is required to solve the equation for t: 

S{t) = 10* +100exp(-0. It) -100 = 100. 

It shows that S (12) = 50.1194 million dollars and it will take about 18.4 months to reach the total sales of $100 
million (the negative solution does not apply). 



This is a Classroom License for instructional use only. 
Research and commercial use is prohibited. 

To get started, select MATLAB Help or Demos from the Help menu. 

» t=1 2; 

»■ 1 0.*t+1 GO.*exp(-.1 *t)-1 OO 
ans = 

50.1 194 

solvef 1 0.*t+1 00.*exp(-.1 *t)-1 00=1 OO') 
ans = 

1 8.414056604369606378466046580125 
-1 1 .461932206205825852370610285214 

Figure 15b MATLAB solution to (2) and (3) in Example 2.8 
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Example 2.9. Consumers’ and producers’ surplus, equilibrium price and equilibrium quantity. These are the 
important concepts in business. 

If (x, p) is a point from the price-demand equation p = D(x ) for a particular product, the consumers’ surplus 
(CS) at a price level p ( p = D(x) ) is defined by 

CS = j [z)(x) - p dx. 

0 

CS represents the total savings to consumers who are willing to pay more than p but are still able to buy the 
product at the price p . 

If (x, p) is a point from the price-supply equation p = S(x ) for a particular product, the producers’ surplus (PS) 
at a price level p ( p = S (x)) is defined by 

RS =j[p-S(x)] dx. 

o 

PS represents the total savings to producers who are willing to supply units at a lower than p but are still able to 
supply units at the price p . 

If p = D(x) and p = iS(x) are the price-demand and price-supply equations, respectively, for a product, and if 
(x, p) is the point of intersection of these two equations, p is called the equilibrium price and X (where 
p = D(x ) = S'(x)) the equilibrium quantity. 

(1) The price-demand equation is given by p = D(x ) —9 — ln(x + 4) . Find the consumers’ surplus at the 
price level p = $2,089 ; 

(2) The price-supply equation is given by p = /S(x) = 51n(x +1) . Find the producers’ surplus at the price 
level p = $25 ; 

(3) Find the equilibrium price and the equilibrium quantity. 

The MATLAB solutions are given in Figure 16 (a, b, c). Figure 16a is for solution of (1), the first command 
tt solve( 4 2.089=9-ln(x+4)’)” is for finding X so that 2.089 = p = D(x) = 9 — ln(x + 4) . It shows that 
X « 999.25 . The second command “int(‘9-ln(x+4)-2.089’,0,999.25)” is for computing the consumers’ surplus 

CS = j[D(x)-p\dx= J [9-ln(x +4)-2.089] ^ = 977.15. 

0 0 

This means that the consumers’ surplus is $977.15. It represents the total savings $977.15 to consumers who are 
willing to pay more than the price $2,089 but are still able to buy the product at the price $2,089. 
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This is a Classroom License for instructional use only. 
Research and commercial use is prohibited. 

To get started, select MATLAB Help or Demos from the Help menu. 

soIve(' 2.089=9-1ogOc-i-4) ') 

ans = 

999.24999082324610099399634909241 

int('9-log(x+4)-2.089 , ,0,999.25) 
ans = 

977.1 51 16826772562149933040020657 

Figure 16a. MATLAB solution to (1) of Example 2.9 


Figure 16b is for solution of (2), the first command “solve(‘25=5*ln(x+l)’)” is for finding X so that 

25 = p = S(x) = 51n(x + 1) . It shows that X = exp(5) — 1 « 147 . The second command “int(‘25- 
5*ln(x+l)’,0,147)” is for computing the producers’ surplus 


X 147 

PS = J[p-S(jc)]<fr = J[25 -51n(jc + l)]dbc = 712. 


06. 


This means that the producers’ surplus is $712.06. It represents the total savings $712.06 to producers who are 
willing to supply at a lower price than $25 but are still able to supply the product at the price $25. 


1 MATLAB 7.3.0 (R2006b) 



File Edit 

Debug Desktop Window Help 



Q 

*<0 CT^n EST 1^3 Current Directory: C:\Program Files\MATLAE\R2006b\work 

- GH 
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Shortcuts 

1 - 

■j" 1 How to Add IS Whet's New 

rtesearcn ana commercial use is Dronioneo 




To get started, select MATLAB Help or Demos from the Help menu. 
» solve('25=5*log(x+1)') 


ans = 
exp (5)-1 

» intC25-5*log(x+1)', 0,1 47) 
ans = 

441 0-1 480*log(2)-740*log(37) 

« 441 0-1 480*log(2)-740*log(37) 
ans = 

71 2.0629 


Figure 16b. MATLAB solution to (2) of Example 2.9 
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Figure 16c is for solution of (3), the first command “solve(‘9-log(x+4)=5*log(x+l)’)” is for finding the equilibrium 
quantity X so that D(x) —9 — ln(x + 4) = S(x) = 51n(x + 1) . It shows that X « 3 . The equilibrium price 
p — 9 — ln(3.09 + 4) « $7.04 per unit. 


1 MATLAB 7.3.0 (R2006b) 


File Edit Debug Desktop Window Help 


O 1 'Jr <K> Hi Gil D ^ ^ ^ Current Directory: C:\ProgramFiles\MATLAB\R2006b\work v [7 

□ * 

Shortcuts \E How to Add \T] What's New 



This is a Classroom License for instructional use only. 
Research and commercial use is prohibited. 

To get started, select MATLAB Help or Demos from the Help menu. 

» solve('9-log(x44)=5*log(x+1 )') 

ans = 


3.0889659101959689709467729053229 


» x=[0:.05:5]; 

» D=9-log(x44); 

» S=5*log(x+1); 

» plot(x,D, , - , ,x,S, , - , ) 

» gtext(' Price-demand eq. D(x)‘) 
» gtext(' Price-supply eq. S(x)‘) 
» gtext('Equilibrium price=7.04') 
» gtext('Equilibrium quantity=3‘) 
» 



Figure 16c. MATLAB solution to (3) of Example 2.9 


2.5. Introduction to multivariate calculus 

Functions of more than one variable are frequently encountered in quantitative business problems. For 
example, if a company produces two types of products (x,y) and each type of products usually has a price-demand 
equation, say, p = p(x , y) and q = q(x, y) , then the revenue from the two types of products is a two-variable 
function 

R(x, y ) = xp(x, y ) + yq(x, y) . 

If the cost function is given as C(x, _v), then the profit function is also a two-variable function 
P(x, y) = xp(x, y) + yq(x, y) - C(x, y). 

If the company wants to maximize its profit based on its available resources, it leads a maximization 
problem on the two-variable profit function P(x , y) . This can be easily extended to the cases of three or more 
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variables if the company produces more than one type of products. Problems related to maximization of multivariate 
profit functions belong to the topic of multivariate calculus. In this section we will illustrate the MATLAB solutions 
to some simple multivariate calculus problems. 

2.5.1. Graphing a two-variable function 

Example 2.10. Plotting the two functions: 

(1 ) Z = f(x,y) = e~ x ~ y2 =exp(-x 2 -y 2 ), -3<x<3, -3<y<3; 

(2) z = f(x, y) = xe~ x ~ y = xexp(-x 2 - y 2 ), — 3 < x < 3, — 3 < y < 3 . 

The MATLAB solutions to (1) and (2) are given in Figure 17a and Figure 17b, respectively. 



This is a Classroom License for instructional use only. 
Research and commercial use is prohibited. 

To get started, select MATLAB Help or Demos from the Help menu. 


» syms x y 

» ezmesh( , exp(-x A 2-y A 2)\ [-3,3,-3,3]) 


i*' „ Figure 1 





File Edit View 

Insert Tools Desktop Window Help 

■■sj 

D £ Si 


<a. ® « 

□ E3 

■ O 


exp(-x 2 -y 2 ) 


11 ^* 



Figure 17a Plot of the two-variable function (1) 
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This is a Classroom License for instructional use only. 
Research and commercial use is prohibited. 

To get started, select MAT LAB Help or Demos from the Help menu. 


» syms x y 

» ezm e s h (V exp (-x A 2-y A 2)' , [-3,3,-3,3]) 



Figure 17b Plot of the two-variable function (2) 


2.5.2. Partial derivatives 


Example 2.11. Find the partial derivatives of the following functions. 

(1) f(x,y) = x 2 + ye x -y 2 -, 

(2) g(x, y, z ) = 2x 2 + xy 2 -3xz 2 ■ 


The MATLAB solutions to (1) and (2) are given in Figure 18a and Figure 18b, respectively. The results in Figure 
18a and Figure 18b show that 


dx 


= 2x+ ye x , 


- = e’-2y 
8y 


8 J- = Ax+y K iy e, ®S=2xy-3« J , 
dx dy 


— = -6 xyz 
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This is a Classroom License for instructional use only. 
Research and commercial use is prohibited. 

To get started, select MAT LAB Help or Demos from the Help menu. 

syms x y 

d i ff 0x^2 +y * exp (x)-y jrtL 2' , x) 
ans = 

2 *x+y * exp (x) 


» d i ff fx A 2 +y * exp (x)-y A 2 1 , y) 

ans = 

exp(x)-2*y 


Figure 18a MATLAB solution to (1) 


MATLAB 7.3.0 (RZOOfib) j 

File Edit Debug Desktop Window 

Help 


X na m « n. | * 

ef E 

Current Directory: C:\Program Files\MATLAB\R2006b\work ■v jT 


Shortcuts 03 How to Add 03 What's New 


i o get started, select i lmd neip or L^emos rrom tne rieip menu. 

>->■ syms x y z 

» diff( , 2*x A 2+x*y A 2-3*x*y*z A 2 , ,x) 
ans = 

4- *x+y A 2- 3 *y *z*2 


» d i ff (' 2 *x A 2 +x*y A 2- 3 *x*y *z*2 1 , y) 
ans = 


2 *x*y- 3 *x*z A 2 


» d i ff (' 2 *x A 2 +x*y A 2- 3 *x*y *z*2 1 , z) 
ans = 

-6*x*y*z 


Figure 18b MATLAB solution to (2) 
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2.5.3. Optimization using Lagrange multipliers 

An optimal strategy for a quantitative business problem is usually subject to finite resources such as the 
finite budget. Therefore, it is a situation that optimization of an objective function (such as the profit function) is 
subject to certain constraints. The Lagrange multiplier method is usually employed to solve an optimization problem 
with constraints. See the following example. 

Example 2.12. A manufacturing firm has budgeted $60,000 per month for labor and materials. If $x thousand is 
spent on labor and $y thousand is spent on materials, and if the monthly output (in units) is given by 
N(x, y ) = 4xy - 8x, 

How should the $60,000 be allocated to labor and materials in order to maximize A? What is maximum A? 

This is a maximization problem subject to the constraint X + y = 60 (in thousands). The objective function using 
the Lagrange multiplier is given by 

fix, y, z) = 4xy-8x + z(x + y -60), 

where z is the Lagrange multiplier. It is well-known that the optimal solution to ( X , y, z) is subject to: 


£-a £-o. 

dx dy dz 

This is related to the computation of partial derivatives. The MATLAB solution is given in Figure 19. 


k MATLAB 

7.3.0 (R2006b) 
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menu. 


» syms x y z 

>*■>*■ di ff('4 x*y-8(x-i-y- 60)',x) 


4*y- 0+z 

» di ff ('<4- *x*y- 8 *x-i-z* (x-i-y- 60)\y) 
ans = 

4*x+z 

>*■ >*■ d i ff('4 *x*y- 8 *x-i-z* (x-i-y- 60)',z) 


x+y-60 


Figure 19a Computation of partial derivatives 
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[x J y J z]=solve( , 4*y-8+z=0 , J '4*x+z=0 , J , x+y-60 = 0') 

x = 

29 


y = 

31 


-116 

Figure 19b (Figure 19a continues) Optimal solution to ( X , y, z) 

Figure 19b shows that the optimal solution to (x, y, z) is x=29, y=31 , and Z = — 116 . Therefore, the maximum 
output is: N = 4xy - 8x = 4 x 29 x 31 - 8 x 29 = 3364. 

In the practical point of view, the above solution to ( X , y,z) will be the best solution to maximizing the 
output function. In the theoretical point of view, it is necessary to check whether the locally optimal solution is 
really the best solution. This is related to the second order partial derivatives. We omit this step in practical business 
applications. 

3. CONCLUDING REMARKS 

MATLAB is a powerful numerical computation package. I t can be used as an interactive mathematical 
shell and/or toolbox. In teaching of calculus-based business mathematics course at University of New Haven in the 
past few years, our experiences have consistently shown us that most freshmen in college of business do welcome 
the opportunity to learn the fundamental business concepts of cost minimization; profit maximization, as well as the 
quantitative techniques of searching for the optimal resource allocation solution in multivariate cases related to the 
real business world problems. Nevertheless, the difficulty is arisen when the concept of limits ; derivative and 
integration rules and formulas in calculus are introduced; the majority of these students often struggle to find the 
correct numerical solution to the problem posed. This lack of algebraic ability and skills, we surmise, is due 
primarily to the previous uneven and the inadequate training and preparation. We hope that by teaching our students 
to utilize a dozen or so simple MATLAB function routines, as a supplementary skill-builder, our students can learn 
calculus-based business mathematics much quicker, better and easier. 

What’s more, It is also our expectation that the simple MATLAB function commands we have 
demonstrated in this paper will enable our students to avoid the trap of boggle down in the memorization of the fine 
detail of mathematical procedure and calculation rules; instead, to focus on the key objective of learning business 
mathematics in the first place — to obtain, with the graphical illustrations on the side, quickly and easily the 
numerical solutions, and thus our student will be empowered to examine the big-picture setting. As a result, these 
students will be more likely to become better decision makers in facing the dynamic and ever-changing business 
challenges. 
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